On the origin of heavy quasiparticles in LiV204 
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An explanation is provided for the heavy quasiparticle excitations in LiV2 04. It differs consider- 
ably from that of other known heavy-fermion systems. Main ingredients of our theory are the cubic 
spinel structure of the material and strong short-range correlations of the d electrons. The large 7 
coefficient is shown to result from excitations of Heisenberg spin i rings and chains. The required 
coupling constant is calculated from LDA+U calculations and is found to be of the right size. Also 
the calculated Sommerfeld- Wilson ratio is reasonably close to the observed one. 

PACS numbers: 71.27.+a, 71.10.-w, 75.40.Cx 



The low-temperature properties of the metal L1V2O4 
show many of the characteristic features of heavy-fermion 
systems like CeAi3 or CeRu2Si2 0]. The coefficient 7 of 
the low-temperature specific heat C = 7T is large, i.e., 
7 = 0.35 — 0.42 J jmolK 2 jlj], depending on samples and 
C(T) shows a broad maximum around 16K. The large 
7 value implies a high density of fermionic low-energy 
excitations. In fact, when C(T) is integrated up to 60K 
it is found that the associated entropy is close to 2R In 2 
where R is the gas constant. This suggests that there 
is roughly one excitation per V ion in the system when 
it is heated up to that temperature, a finding inconceiv- 
able for a conventional, partially filled conduction band 
of d-electrons. In fact, the 7 value obtained from ab ini- 
tio band structure calculations |^|-[?J turns out to be too 
small by a factor of 25. In accordance with the observed 
large 7 value it is found that also the spin susceptibility 
Xsp(T) is similarly enhanced ||. It is nearly T indepen- 
dent below T ks 30K with a shallow broad maximum 
around 16K and a sample dependent increase below 5K . 
The Sommerfeld- Wilson ratio Rw = ^k^Xsp/i^P-lel) is 
1.7 with considerable error bars because of the remaining 
low temperature variation of Xsp(T). A ratio of that size 
shows that the enhancements in 7 and Xsp nearly can- 
cel each other as expected for a metal with an enhanced 
density of states N*(0) at the Fermi energy. Remem- 
ber that 7 ~ N*(0) and Xs P ~ iV*(0) in that case. No 
(anti)ferromagnetic order was found above 0.02K and no 
evidence for spin-glass behavior. There are also no indi- 
cations of charge ordering, i.e., all V sites are equivalent. 
The low-temperature resistivity is p(T) = po + AT 2 with 
a large coefficient A — lySlcmj K 2 like in other heavy- 
fermion systems. L1V2O4 is therefore a rare case of a 
system with heavy quasiparticle excitations which do not 
involve / electrons Q. 

The physical origin of the low-energy scale associated 
with the characteristic heavy-fermion excitations has re- 
mained unclear. Presently we know of a few physical 



mechanisms resulting in heavy quasiparticles which all 
seem inapplicable here. 

One well known source is the Kondo effect It ap- 
plies to heavy-fermion systems like CeAl3 and requires 
ions with quasi-localized / electrons of nearly integer 
number. In fact 7 is the larger the closer the /-count 
rif is to an integer, e.g., to n/ = 1 in the case of Ce ions. 
But in L1V2O4 the valency of the V ions is +3.5 and 
therefore the Kondo effect is excluded here. A second 
mechanism is realized in Nd2- x Ce :c Cu04 (x « 0.1 — 0.2) 
and is essentially based on the Zeeman effect ||. It ap- 
plies to systems with strongly correlated conduction elec- 
trons which couple weakly to ions with localized spins, 
i.e., Nd ions in the above case. This situation can be 
also excluded for LiV2 04. The same holds true for heavy 
fermions appaering very close to a quantum critical point, 
a model which has been suggested for YMn2 [|l0| . There 
remains the case of heavy quasiparticles caused by charge 
ordering like in Yb4As3. Here the characteristic fea- 
ture is the formation of well separated spin i Heisenberg 
chains in a three-dimensional lattice due to charge order 
p],^2| . This route seems promissing except that there 
is no charge ordering in L1V2O4. In that material all V 
sites remain equivalent down to the lowest temperatures. 
Nevertheless, it is this example which leads us to suggest 
the physical model discussed in the following. 

The aim of this communication is to discuss a new sce- 
nario for the origin of a low-energy scale. It uses the 
fact that L1V2O4 has a cubic spinel structure down to 
the lowest temperatures measured P,fl3"f , i.e., the sys- 
tem seemingly does not charge order. But it has ob- 
viously strong short-ranged electron correlations. The 
spinel structure enables us to account for them in a par- 
ticular way which suggests directly a plausible explana- 
tion of the large number of fermionic low-energy excita- 
tions as evidenced by the large 7 coefficient. 

In the spinel structure the V ions form corner-sharing 
tetrahedra. Their sublattice can be viewed as a fee lat- 
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tice with every second site removed. A V ion has six 
nearest neighbors of its kind and the V ions are the cen- 
ters of slightly distorted edge sharing octahedra of oxy- 
gen ions. All V ions are crystallographically equivalent. 
Spinel structures or pyroclore lattices give raise to frus- 
trated antiferromagnetic interactions. Resulting forms of 
a spin- liquid ground-state have been frequently discussed 
see, e.g., 0||. 

Band structure calculations employing the LDA+J7 
computational method show that the Fermi energy is 
placed in vanadium t2 g bands which are well separated 
from the oxygen p-bands. The average d-electron num- 
ber per V ion is 1.5. For those electrons we assume the 
following simplified model Hamiltonian to apply 

H = ^2 t "( a ro(») fl w(i)+' l ' c ') 
(ij),va 

+ U y j n v -\ [i)n v i (i) + U }^ n„(z)n M (i) 
+ V Q ^2n(i)n(j) + J H ^2s 1 (i)s 2 (i) 

(ij) ' 

+ £jyS(i)S(j) . (1) 

(«> 

The a+ .(i)(a J/cr (i)) create (destroy) a <i-electron in orbital 
v with spin a at site i. For simplicity we shall use two d- 
orbitals only, one pointing along the octahedron axis and 
the other pointing towards the four oxygen ions within 
the plane. The first term describes the kinetic energy and 
only transitions between neighboring sites (ij) are taken 
into account. The second and third term account for 
the on-site Coulomb repulsion of <i-clcctrons in the same 
and in different orbitals. Here = ^2 aJ .(i)a y[T (i) 

a 

is the number operator. The fourth term describes the 
Coulomb repulsion of electrons situated on neighboring 
sites, the fifth term is due to Hund's rule coupling and 
the last terms results from superexchange via the oxygen 
sites. The total spins S^tj) of the neighboring sites i and 
j depend on the number of d electrons at those sites. For 
a d 1 configuration S — | while for a d 2 site S = or 
1 etc. We assume that U and Vq are much larger than 
the t v while Jij is expected to be of similar size. A large 
value of U implies that only d 1 and d 2 configurations need 
consideration. Moreover, a large value of Vq requires that 
on each tetrahedron two sites are in a d}(y A+ ) and two 
are in a d 2 (V 3+ ) configurations (tetrahedron rule). This 
has been first pointed out by Verwey |t7| and was dis- 
cussed in detail by Anderson jio]]. Note, that because of 
an average d occupancy d 2 — d 2 and d 1 ~ d 1 configu- 
rations repel each other while d 1 — d 2 configurations are 
attractive. LiV204 is a metal and therefore the above 
tetrahedron rule cannot be strictly obeyed. Otherwise a 
single transition d 2 =^ d 1 would imply a rearrangement 
of configurations in nearly all tetrahedra. However, we 
are interested here in the spin excitations only, because of 



the large entropy associated with the low-energy excita- 
tions. Therefore the heavy quasiparticle excitations must 
involve predominantly spin degrees of freedom. Thus for 
the present purpose we may assume that the tetrahedron 
rule is strictly obeyed. In that case the ground state | Vo) 
of the system is of the form 

l<M=5>n|<M (2) 

n 

where the \<j) n ) are all configurations which obey the 
tetrahedron rule. In each of the \cf> n ) the V A+ ions form 
chains and rings. They are surrounded by chains and 
rings of V 3+ ions. An example is shown in Fig. 1. One 
may check that the smallest rings possible involve six 
d l {d 2 ) ions. Ring sizes increase in units of four, i.e., they 
involve 10,14 etc. V ions of one kind. The form of the 
a n is irrelevant here and will be the subject of a separate 
investigation. 

Consider the d 2 ions. Without Hund's coupling the 
ratio of S = 1 triplet and S = singlet configurations 
would be approximately 3 : 1. Finite Hund's coupling 
enlarges it further. Therefore we deal with 5=1 chains 
(and rings) interrupted occasionally by sites with S = 0. 




FIG. 1. Example of a configuration \4>n) obeying the tetra- 
hedron rule. Solid dots: spin i sites; empty dots: spin 1 sites. 
They form chains and rings each. Notice a ring consisting of 
six spin \ sites in the upper right corner. 



Let us assume next that the coupling between S — 1 
and S — sites would vanish. In that case the exci- 
tations of the system are given by those of the S = \ 
chains and rings and of the S — 1 chains and rings, the 
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latter being disrupted by 5 — configurations. At low 
temperatures only the excitations of the 5 = | Heisen- 
berg chains are important because their spectrum goes 
to zero in the long wave length limit. The excitations of 
the 5 = 1 Heisenberg chains and rings have a gap in the 
excitation spectrum and play no role at low T. The same 
holds true for 5 = \ rings which are sufficiently small. 
The relevant excitations are of the form 



(3) 



where A\ generates a spin wave of wavelength A = =£■ 
i chains and large rings contained in \(j) n ). 



in the 5 = 1 



The 5 



chains and large rings have a much larger 



statistical weight than the 5 — % small rings and there- 
fore we may assume that half of the N sites are involved 
in the low energy excitations. The corresponding low- 
temperature specific heat is C = 7T with a 7 coefficient 
given by § 



7=o 



2k B R 



3 J 3 



(4) 



where J3 = when S(i) = S(j) = ^. Note that 
only one V ion per formula unit is contributing to 
the low-energy excitations and hence to 7. The low- 
temperature spin susceptibility is that of Heisenberg 
chains and given by [p| Xsp — ^^IgR/i^ 2 J3) which re- 
sults in a Sommcrfcld- Wilson ratio of Rw = 2. 

In the following we want to provide arguments which 
suggests that the above considerations still apply, when 
the coupling Jy between sites with spin 5 = \ and 5 = 1 
is included. For that purpose we must first estimate the 
relative sizes of Jij when spins \ and spins 1 are cou- 
pled among themselves, i.e., J3 and J\ and when a spin 
\ couples with a spin 1, i.e., J2- With this aim we have 
performed LDA+J7 band structure calculations in order 
to estimate their relative sizes. The calculations were 
performed under the following assumptions. First of all, 
we assumed that the V ions with d 1 and d configurations 
form two different families of ordered chains which alter- 
nate in z direction. Both families lie in the xy plane but 
they differ by a rotation of tt/2. With this arrangement 
each V tetrahedron is occupied by an equal number of d 1 
and d 2 ions, i.e., the tetrahedron rule is strictly obeyed. 
It should be noticed that no charge ordering of V ions 
has been observed experimentally. Nevertheless the cho- 
sen configuration serves its purpose. We have found that 
for such an arrangement of V ions a self-consistent so- 
lution can only be obtained for values of the screened 
on-site Coulomb integral U larger than 3.8 eV, where as 
the use of smaller values of U leads to a solution with the 
average d 15 count for each V ion. The value of U = 4 
eV, that was used in the present work to calculate , is 
larger than the value of about 3 eV previously estimated 
and the experimental one which is close to 2 eV P,pO|. 



However, we believe that the calculated Jij can be used 
to estimate their relative strengths. In order to calculate 
the coupling constants we started from a self-consistent 
solution for a collinear magnetic state with an AF or- 
dering of d 1 and d 2 ions along the corresponding chains. 
In this case each d 1 ion has two d 2 nearest neighbours 
with the same and two with the opposite orientation of 
the spin direction, i.e., d 1 — d 2 magnetic interactions are 
frustrated. Then, the band energies of non-collinear mag- 
netic structures with an angle between either the d 1 spins 
or the d 2 spins varying from 180° to 0° were calculated, 
and Jij's were determined by comparing total energies 
with those of the corresponding Heisenberg-likc model 
with nearest neighbour interaction. (For a detailed de- 
scription of the procedure see Ref. Q and references 
therein.) As a result we obtained values of 3 meV and 24 
meV for the coupling constants J3 and J\, respectively. 
The coupling of a spin i with a spin 1 is intermediate. 




FIG. 2. Two tetrahedra of a spinel structure. Sites 
v = 1, 2, 3 have spin S(y) = 1 while the remaining ones have 
spin i 

Consider first two corner-sharing tetrahedra (see Fig. 
2) with spins S(z/). Let sites 1, 2, 3 have spin 5=1 
while the remaining sites have spin 5 = i . We introduce 

Ej = S(l) + S(2), E 2 = S(2) + S(3), u x = S(4) + S(5), 
&2 = S(6) + S(7) in terms of which the interaction Hamil- 
tonian Yluj) is written 
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£l<72 + E2CT1 



(5) 



The ground state has = = o\ = a% = 0. 
When Ji J2, J3 the low-energy excitations are given 
by changing o\ or er 2 from to 1, i.e., they are within the 
spin i subsystem. They are unaffected by J2, i.e., the 
coupling between spin ^ and spin 1 sites. The argument 
builds on the fact that the coupling J\ between the spin 
1 sites is much stronger than the one between the spin 
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\ sites, i.e., J3. It should hold also when instead of two 
tetrahedra the whole lattice is considered. The excita- 
tions of spin I chains and large rings are effected by the 
5 = 1 subsystem only via spin flips of the latter. How- 
ever, the excitations of the S = 1 subsystem have a gap 
and moreover from the LDA+C/ calculations it follows 
that Ji ~ 8J3. Therefore the low- lying excitations take 
place within the S = \ subsystem only. Consequently 
7 is given by that of uncoupled S = h chains (or large 
rings), i.e., by (|4|). Indeed, this reasoning provides for a 
simple explanation of the heavy quasiparticle excitations 
found in LiV2 04. In order to explain the experimentally 
observed value of 7 = QA2J/mol ■ K 2 a coupling con- 
stant Js/ks — 13. 3K is required. This value agrees up 
to a factor of 2.6 with the value J| alc //cs = 35K found 
from the LDA+[7 calculations. Furthermore, the ob- 
served Sommerfeld- Wilson ratio of 1.7 is relatively close 
to the theoretical value of Rw — 2 found here. Some 
comments are in order at that point. We cannot expect 
that a LDA+l/ calculations give precise values of spin 
coupling constants, since those physical quantities are 
rather sensitive to a precise description of correlations 
(see, e.g., p2| ) which would require much more sophisti- 
cated computations. Also 7 as given above may well con- 
tain contributions from phonons, i.e., 7 = 70(1 + A) where 
A describes the contribution from the electron-phonon in- 
teractions. For example, a value of A = 1 would double 
the value of J3 needed for explaining 70 to J3 27K. 
But it would also double the experimental value of Rw 
to 3.4 since A does not contribute to Xsp- The calculated 
value of Rw = 2 is a result for the spin ^ subsystem. 
Therefore any contribution of the spin 1 subsystem |lTj| ] 
to the measured values of \s P must be subtracted before 
a comparism with experiments is made. Even a small 
contribution of 15 % reduces Rw from 3.4 to 2.7. 

Experiments show no maximum in j(T) — C(T)/T 
in distinction to theor y |18| and a broad one in x sp (T) 
around T = l&K [p] j2^]25 ]. Theory predicts such a max- 
imum at 0.6 J3 which would imply a value of J3 ~ 27 'K . 

Although the agreement between theory and experi- 
ment is not perfect, which comes to no surprise in view 
of the simplifications made, we believe to have found a 
simple physical picture for the hitherto unexplained un- 
usual low-temperature thermodynamic properties of the 
material. 

It is based on the proposition that the system is a 
superposition of configurations with spaghetti like spin 
i chains and rings. They are well separated from each 
other since the spin 1 chains and rings surrounding them 
have a large coupling constant J\ and a gap in the excita- 
tion spectrum (Haldane gap). Therefore, the low-energy 
excitations take place in spin h Heisenberg chains only 
similarly as in YhiAs3. Because of the near separation 
of spin and charge degrees of freedom LiV204 should not 
be a Fermi liquid in the classical sense. Rather it should 
show features of a Luttinger liquid because of the pres- 



ence of an intrinsing chain structure even in the cubic 
lattice. It would be most interesting to determine the ef- 
fective masses at the Fermi surface, e.g., by de Haas- van 
Alphen experiments. We speculate that they should be 
much smaller than suggested by the large value of 7. As 
pointed out above the tetrahedron rule is not expected 
to hold rigorously. Therefore, the influence of small devi- 
ations from it remains an interesting topic which will be 
investigated separately. The same holds true for the ef- 
fect of impurities and moreover for doping. Furthermore, 
it remains to be seen whether similar ideas apply also to 
LiTi 2 04 which is a superconductor at low temperatures 
or to pyrochlore lattice systems. 
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